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Inspired by a recent experimental observation of the zero-bias tunneling conductance in 
superconductor-semiconductor nanowire devices [T], we consider here transport properties of the 
junctions consisting of a nanowire (Luttinger liquid) coupled to a topological superconductor char- 
acterized by the presence of Majorana zero-energy end states. The presence of the Majorana modes 
leads to a quantization of the zero-bias tunneling conductance at zero temperature. In order to un- 
derstand this phenomenon, we have developed a framework, based on real-time Keldysh technique, 
which allows one to compute tunneling conductance at finite temperature and voltage in a realistic 
experimental setup. Our approach allows one to understand this transport phenomenon from a more 
general perspective by including the effect of interactions in the nanowire, which sometimes results 
in a drastic departure from the non-interacting predictions. Thus, our results provide a key insight 
for the tunneling experiments aiming at detecting Majorana particles in one-dimensional nanowire 
devices. 

PACS numbers: 73.21. Hb, 71.10.Pm, 74.78.Fk 



I. INTRODUCTION 

The search for topological phases in nature 
has been an active and exciting pursuit. Re- 
cent experimental progress [iHS] aimed at de- 
tecting Majorana zero-energy modes (Major anas) 
in one-dimensional superconductor /semiconductor 
nanowire devices has excited the whole physics com- 
munity . The reason for this excitement is two- 
fold. First of all, Majorana modes are predicted 
to obey non-Abelian braiding statistics [9HT3] and, 
thus, their discovery is of fundamental importance 
for all of physics. Second, this property of Majo- 
ranas (or Ising anyons) is at the heart of topolog- 
ical quantum computing schemes [HI [15]. Indeed, 
the Majorana modes localized at the opposite ends 
of the nanowire provide a one-dimensional topologi- 
cally protected two- level system [16 . The network of 
Majorana wires can be used for topological quantum 
information processing p!7''^20^. The topological ap- 
proach to quantum computing is very promising be- 
cause it provides a way to protect information from 
local errors by encoding it in non-local degrees of 
freedom of topologically ordered systems [HI |15] . 

It has been predicted that Majorana zero-energy 
modes can emerge quite naturally in spinless p- 
wave superconductors [11] [16]. Experimental real- 
ization of a spinless p-wave superconducting state, 
which is robust against various perturbations inher- 
ently present in any solid-state system (e.g. dis- 
order), is a non-trivial task. One attractive class 



of experimental proposals involves a conductor with 
strong spin-orbit coupling allowing one to enslave 
the spin degree of freedom and an s-wave supercon- 
ductor [2TH26] . The experimental search of Majo- 
ranas in solid-state systems has been focused so far 
on one-dimensional superconductor/semiconductor 
nanowire devices ^3 HI] because of the simplicity 
and robustness of these systems. (We refer a reader 
to the following review articles [27l [28i for a de- 
tailed discussion and comparison of different exper- 
imental proposals.) The presence or absence of Ma- 
jorana zero-energy modes in one-dimensional topo- 
logical superconductors can be detected by measur- 
ing tunneling local density of states at the ends of 
the wire. Shortly after the Majorana wire propos- 
als [25] [26], there have been several experimental 
reports [H [SHS] of zero-bias tunneling conductance 
in superconductor/semiconductor nanowire devices, 
which can be interpreted as a signature of Majorana 
zero-energy modes. This interpretation of the ob- 
served zero-bias tunneling conductance relies on the- 
oretical predictions based on non-interacting models 
for one-dimensional semiconductor-superconductor 
heterostructures [29ti38 . 

In this paper we revisit transport properties of 
a TP SC/LL junction from a much more general 
perspective using a combination of renormalization 
group (RG) and real-time Keldysh techniques. Our 
approach allows one to take into account electron- 
electron interactions, which are particularly impor- 
tant in one-dimensional systems [39], and com- 
pute transport properties of the superconductor- 



semiconductor nanowire devices as a function of var- 
ious parameters such as temperature, voltage, and 
length of the nanowire. The system of interest in- 
volves a topological superconductor coupled to a 
Luttinger liquid representing the uncovered portion 
of a nanowire, see Fig. [l^ for the schematic plot of 
the experimental setup [1 . Thus, the simplest the- 
oretical model for the tunneling experiment [1 con- 
sists of a topological superconductor coupled to a 
one-dimensional nanowire of length L/2. The other 
end of the nanowire is connected to a semi-infinite 
normal-metal lead which is a part of the circuit for 
measuring tunneling conductance. Here the metallic 
lead is described as a three-dimensional Fermi liquid 
where interaction effects between quasiparticles can 
be neglected. For our purposes, it can be represented 
as a non-interacting Luttinger liquid with Kl = 1. 
Therefore, in order to study transport properties of 
such a system, we introduce an inhomogeneous LL 
interaction parameter K{x) which adiabatically in- 
terpolates between K{x) = Ky^ and K{x) = Kl in 
the nanowire and leads, respectively. The schematic 
plot of the experimental setup is shown in Fig. ^p. 

In the absence of the lead, the properties of TP 
SC/LL junctions have been studied in Ref. [40] 
where it was shown that the quantization of the tun- 
neling conductance in units of 2e^ /h is robust even 
in the presence of electron-electron interactions in 
the wire. In particular, for weaker electron-electron 
interactions {Ky^ > 1/2), the system fiows to the so- 
called perfect Andreev reflection fixed point where 
the zero-bias tunneling conductance is equal to the 
conductance quantum at zero temperature. We note 
here that in the case of a non-topological supercon- 
ductor coupled to a Luttinger liquid [4TH43] , the per- 
fect Andreev reflection fixed point is unstable for any 
repulsive interparticle interactions in the wire, i.e. 
for Kyj < 1. Thus, there is a qualitative difference 
between TP SC- and NTP SC-LL junctions. 

Building on the previous results [40 , we study 
here the transport properties of a TP SC/LL junc- 
tion in the presence of the leads. Specifically, we 
derive the Keldysh action for the boundary theory 
describing such a junction and then compute its 
transport properties. We first show unambiguously 
that tunneling conductance is quantized in units of 
2e^//i, as expected from the previous studies [44- 
[54] , Next, we compute temperature- and voltage- 
dependent corrections to the tunneling conductance 
and discuss their scaling in the presence of the leads. 
We show that in some instances the presence of 
interactions in the nanowire results in a qualita- 
tively different dependence of the tunneling conduc- 
tance on physical parameters compared to the non- 
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FIG. 1. a) Schematic picture of the semicon- 
ductor/superconductor hybrid device in Delft experi- 
ment fP. b) Schematic picture of the theoretical model 
for the corresponding device. 



interacting predictions. Thus, our results provide 
important information for the Major ana tunneling 
experiments, particularly since the maximum zero- 
bias peak value is smaller than 2e^ /h by an order 
of magnitude in current tunneling experiments [T, ^ 
5 . The nature of the reduced tunneling conductance 
from the quantized value 2e^ /h at finite temperature 
and voltage is an open question right now. 

The theoretical framework developed in this pa- 
per has broader impacts and can be applied to 
other quantum impurity problems, see, e.g. , lec- 
ture notes [55] and references therein. We note 
that the presence of the leads breaks integrability 
of the boundary sine-Gordon model, and the ex- 
act solution [56 is not available in this case. Thus, 
perturbative calculations such as the one developed 
here are particularly useful. Furthermore, the real- 
time Keldysh approach is suitable for studies of 
time-dependent phenomena. Indeed, it should be 
straightforward to compute noise correlation func- 
tions using this method. 

The paper is organized as follows. In Sec. |TTj 
we briefiy review previous results for the tunneling 
conductance in one-dimensional Majorana wires. In 



Sec. IIL we introduce Keldysh formalism and com- 
pute transport properties of a TP SC/LL junction 
first for an infinite LL and then for a finite LL cou- 
pled to a metallic lead as shown on Fig. [1J3. In 
Sec. IV, we summarize our main findings and pro- 
vide a qualitative explanation of the temperature 
and voltage dependence of the conductance using a 
two-step RG procedure. The technical details are 
relegated to the Appendices. 



2 



II. DISCUSSION OF THE PREVIOUS 
RESULTS 

A. Zero-bias anomaly and tunneling 
conductance in the non-interacting limit 

In this section, we briefly review previous results 
for transport calculations in the non-interacting 
model of TP SC/LL junctions [29ti38]. The pres- 
ence of the lead is not important here since Ky^ = 
Kl = 1, and we assume that L ^ oo. 

Without loss of generality, we focus here on the 
transport properties of a spinless nanowire coupled 
to a spinless ID superconductor. For simplicity, we 
consider below a single-band nanowire but our re- 
sults should also be applicable to multi-band super- 
conductor/semiconductor heterostructures [57-59 . 
Throughout this paper, we will assume that the su- 
perconductor is fully gapped and sufficiently long 
so that the splitting between Majorana end states 
can be neglected. Therefore, the LL is effectively 
coupled only to 71 which is the only relevant low- 
energy superconducting degree of freedom in the 
problem. Taking into account the aforementioned 
assumptions, the effective model for the topological 
superconductor-LL junction becomes {h = 1) 



H = H, 



(1) 
(2) 



where v is the Fermi velocity and i^l^^j^ represents 
right/left-moving excitations near the Fermi energy. 
It is convenient to perform an unfolding procedure 
and rewrite Hq in terms of a single fermion field ip{x) 
defined over all X as follows: 



^pnix), x>0 



In terms of '^(x), Hq becomes simply 



dx (^—ivip^ dxi^) 



(3) 



(4) 



The Hamiltonian SH in Eq.Q represents bound- 
ary terms. As shown in Ref. |40 , the most impor- 
tant processes from the RG perspective are tunnel- 
ing to Majorana mode and electron backscattering 
from the boundary defined by the amplitudes Am 
and A AT, respectively. The boundary Hamiltonian is 
given by 



6H 



-f 

J — ( 



^)^2\n^^^]8{x). (5) 



The tunneling current for this setup can be calcu- 
lated using Landauer formalism. The subgap trans- 
port at low temperature and voltage T, eV <C A 
is determined by the Andreev reflection probability 
T{E) [271 EH]: 



2e 
~h 



dE{fF{E-eV)-ME))T{E), (6) 



where Jf^E) is the Fermi distribution function. A 
straightforward calculation of the Andreev reflection 
probability T{E) yields 



T{E) = 



r2 + ^2 ^ 



(7) 



where F = v{\m /vpf /{l^ {\n /vpf), see Ref. gO] 
for details. The differential tunneling conductance 
G = dl/dV can be calculated in various limits using 
Eqs.(|6| and At zero temperature T = 0, one 
arrives at 



G 



{evy 



2e^/h 



(8) 



1 



ieVf 

r2 J 



eF<cr. 



At finite temperature T ^ eV^ the conductance G 
is given by 



G 



2e^/h 



^ T ^ F 



i-^,T<r. 



(9) 



One can notice that the scaling of the conduc- 
tance with temperature and voltage in the limit 
max{eF, T} F is different, which is rather puz- 
zling since the scaling theory predicts the same 
power-law dependance for temperature and volt- 
age [40]. We will show below how to resolve this 
contradiction. It turns out that the noninteracting 
result is non-generic and as soon as we include in- 
teractions (i.e. 1), the scalings of the conduc- 
tance with T and eV become the same. 



B. The phase diagram of a spinless TP SC/LL 
junction 

In this section, we discuss how interparticle inter- 
actions in the nanowire affect transport properties 
of the TP SC/LL junction. Interaction effects in 
the topological superconductor itself have been pre- 
viously discussed in Refs. [6QH64] . The main conclu- 
sion of these publications is that weak interparticle 
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FIG. 2. Phase diagram for TP SC/LL junction as a 
function of LL liquid parameter in the wire K. 

interactions do not destroy the topological supercon- 
ducting state but rather lead to the renormalization 
of some physical parameters, such as induced coher- 
ence length. Therefore, the simplified description of 
a Majorana wire in terms of two delocalized zero- 
energy modes captures the basic physics of an inter- 
acting topological superconductor. Consequently, in 
the remainder of the paper, we focus on the effect of 
interactions in the LL which can significantly affect 
transport properties of the TP SC/LL junction. We 
first briefly review the phase diagram for topological 
SC/LL junctions obtained in Ref. [40]. Throughout 
this section, we assume for simplicity that the length 
of the wire ^ oo in which case one can use the 
RG approach to obtain the phase diagram. A major 
virtue of the RG approach is that it allows one to ex- 
tract universal transport signatures in the presence 
of interactions by computing the scaling dimension 
of various perturbations. Thus, one can make rather 
general statements about transport properties of TP 
SC/LL junctions. 

It has been shown in Ref. [40 that there are two 
fixed points to which the system renormahzes de- 
pending on the LL parameter see Fig|2l These 
two fixed points correspond to perfect normal and 
perfect Andreev reflection and represent two different 
boundary conditions at the junction (i.e. at x = 0). 
The former corresponds to tljR{x = 0) = V^l(^ = 
0) which, according to the bosonization convention 
i^R/L{x) ^ e^iGix)±ci){x)) results in the pinning of 
the 0-field at the boundary, i.e. (j){x = 0) = Omodyr. 
The corresponding imaginary-time boundary action 
at the normal-reflection fixed point is defined in 
terms of the dual fluctuating field 6: 

Sn = ^I JklloMI^ (10) 

Here K is Luttinger liquid parameter characterizing 



the strength of the interactions and Q{r) = 0{r^x = 
0). At the perfect normal refiection fixed point, the 
current through the system is zero since electrons 
are completely backscattered from the junction. 

In the opposite limit (i.e. at the perfect Andreev 
reflection fixed point), the field O is pinned at the 
boundary and the corresponding low-energy action 
is given by 




where ^{r) = <p{r^x = 0). In this limit, the electric 
current is finite at zero temperature and is carried by 
electrons undergoing resonant Andreev reflection. 

We now discuss the stability of different fixed 
points in the presence of interparticle interactions. 
Following the RG analysis [40j, one finds that the 
perfect Andreev refiection fixed point is stable for 
K > 1/2. In this domain, the leading relevant opera- 
tor corresponds to the Major ana-coupling boundary 
term 

Sm = 2Am J dra^'cosQ. (12) 

Here cr^ is the Pauli matrix associated with the two- 
fold ground state degeneracy of the topological su- 
perconductor. As shown in Ref. [40 , cFx is a con- 
served quantity, i.e. it has no imaginary time dy- 
namics and for our purposes can be ignored. 

To understand the stability of the perfect normal 
and perfect Andreev reflection fixed point, we now 
analyze the effect of various perturbations. Start- 
ing at the normal-reflection fixed point, one finds 
that the leading relevant perturbation corresponds 
to tunneling onto the Majorana zero-energy mode. 
The RG flow of the coupling constant Am reads 

Thus, Majorana tunneling destabilizes the perfect 
normal reflection fixed point for > 1/2 and drives 
the system towards the perfect Andreev reflection 
fixed point. In this domain of the LL parameter 
the tunneling conductance approaches the famil- 
iar quantized value G = 2e^ /h at zero temperature 
and zero bias even in the presence of the arbitrary 
symmetry-allowed boundary couplings. 

In the opposite limit K < 1/2, the perfect An- 
dreev reflection fixed point becomes unstable due 
to the normal reflection at the junction. This pro- 
cess corresponds to the most relevant operator for 
K < 1/2. The corresponding imaginary-time action 



4 



at the boundary reads 

Sn = -\n J dTC0s(2$). 
Indeed, the RG flow for A at reads 



dX 



N 



dl 



{1-2K) \n. 



(14) 



(15) 



As expected, one finds that the tunnehng conduc- 
tance in this domain scales to zero at zero tempera- 
ture and voltage. The results concerning the stabil- 
ity of different fixed points are summarized in Figj2j 



III. TUNNELING CONDUCTANCE 
CALCULATION USING REAL-TIME 
KELDYSH FORMALISM 

To make a connection with tunneling experiments, 
it is important to take into account the corrections 
to the conductance due to finite temperature and 
voltage, finite length of the nanowire L/2, as well as 
the presence of the leads, see Fig. [T] Henceforth, we 
focus on the experimentally relevant regime K > ^ 
where the perfect Andreev reflection fixed point is 
stable. Throughout the paper we will use Keldysh 
formalism (see Appendix [A| for details and our con- 
ventions) which, as we argue below, is particularly 
suitable for the problem at hand. First of all, it al- 
lows one to maintain the analytical structure of the 
perturbation theory at every step of the calculation 
and avoids the subtleties with analytical continua- 
tion appearing in Matsubara formalism. Second, one 
can distinguish unambiguously the frequency, tem- 
perature, and voltage dependence of the response 
functions, which enables one to go beyond the scal- 
ing theory and obtain the conductance as a function 
of an arbitrary ratio between voltage and tempera- 
ture. 

In the spirit of the Kane-Fisher approach [45l [65] 
[66] . we first integrate out bulk LL modes and de- 
rive a boundary sine-Gordon model. We then as- 
sume that the system renormalizes to a particular 
critical point and study the effect of the leading rel- 
evant/irrelevant perturbations in the vicinity of the 
critical point. This allows one to understand the 
transport properties as a function of voltage and 
temperature as well as other physical parameters. 



A. Conductance near perfect normal reflection 
fixed point 

The relevant model for the TP SC/LL near the 
unstable normal reflection fixed point reads 



DGexp 



K 

"2^ 



duo 
2^ 



uj\\e{uj)\' 



(16) 



-^^M J dr COS ^0(r 



a(r) 



Here a(r) represents the gauge field (i.e. eV{t) = 
in real time). The current through the junction 
is given by /(r) = 6 In Z/Sa{r)[65 . 

We proceed with the Keldysh calculation of the 
tunneling conductance using the boundary model 
defined above and introducing Keldysh contour, see 
Appendix |A] After Keldysh rotation of the fields 
©d/g = 0+ i 0- the corresponding real-time 
action is given by 



— 2zAm 



-SolOci^eq] 



(17) 



20 ci (t) + aci it)\ ^ 26, (t) + a, {t) 



rjdtsml^ 4 j-'y 4 

where 5*0 [6c; , 6,] is defined as 

So[Qci,&c,] = J ^(eci,eq)^GeiH(eci,eq)_^ 

(18) 

with the Green's function G'^^{u j) b eing the 2x2 
matrix in Keldysh space, see Eq.(A5). The current 



through the junction can be obtained using the fol- 
lowing relation 



{Ici{t)) = 2ie 



_5Z_ 
Saq{t) 



|ag(t)^0- 



(19) 



By explicitly taking the functional derivative and 
computing the response function to the lowest non- 
vanishing order in Am, one finds that current is given 

by 



(Iclit)) 



dt'Ce{t-t') 



(20) 



where the connected correlation function Ce(^ — t') 
reads 



5 



Ceit-t') = {sin 
1 

= - exp 



2eci{t)^aci{t) 



sm 



2eci{t')^aci{f) 



cos 



sm 



(21) 



sm 



sinh 



sinh 



with J0^^^ = 0(t)±0(t'). Here the average is taken 
with respect to So[&ch ©g]- 

We win assume throughout this section that the 
thermal coherence length Lt = v/T is shorter than 
the length of the LL {Lt ^ L) in which case one 
can ignore effect of the leads, see Sec jlll C| for de- 
tails. With these approximations, the correlation 
functions {5Gli) and {Sejy^ Sej"^^) can be calculated 
analytically, see Eqs.( |B8| ). Next, given that in the 
dc limit aci{t) — aci{t') is antisymmetric in t — t\ 



the expression (21) can be simplified by dropping 
the antisymmetric part of the integrand. The cor- 
responding expression for the tunneling current is 
given by 



/CO 
dt^ exp 
-oo 



(22) 



. ' aci{t)-aci{t') 
X sm — sm 



2K 



sgn(t-t') 



ttT \ A J 



'f 

Jo 



dx ^^sin( 



/ TT 



sinh^(x) ^2K J 
(23) 



Here A is an ultraviolet cutoff. As expected, the 
expression for current in the dc limit is independent 
of time. After evaluating the integral (23), the dc 
current through the junction is given by 



^47rT 



7rT\ ^ 1 ^ 



1 

K 



sm 




(24) 



where V 



eV 
47rT 



and r(x) is the gamma function. 
Eq.(|24[) is the main result of this section as it allows 



one to study current as a function of an arbitrary 
ratio of two energy scales eV and T. The I — V curve 
non-trivially depends on the interaction parameter 
K and in general (i.e. for 7^ 1) is non-Ohmic. 

One can now obtain the tunneling conductance 
G = dl/dV in different limits. We first discuss the 



eV <C T case. After straightforward manipulations, 
the dc conductance reads 



G 



M 



h A2 



7tT\ ^' 

\2K) 



(25) 



1 
2K 



As one can see, the conductance G is essentially de- 
termined by the scaling dimension of Am which fol- 
lows from the RG flow. The prefactor fi{K) is fi- 
nite in the non-interacting limit fi{K = 1) ^ ir^/S. 
Thus, the result (25) is consistent with the non- 
interacting limit. 

The non-linear tunneling conductance in the limit 
> T reads 



G = 



f2{K) = 



h \A 
_4^shi^ 
2* 



r 2 



f2{K) 



1 

K 



(26) 
(27) 



Similarly to the previous case, the power-law depen- 
dence of G (X Vk~'^ is consistent with the scaling 
theory. However, the pre-factor /2(^) ^ 7r^(l — 
K) vanishes in the non-interacting limit K = 1 
and, thus, the conductance is zero at this order 
of perturbation theory in Am (up to corrections 
0(exp[— eF/T])). This is consistent with the non- 
interacting results discussed in the previous section 
where we showed that G oc cf. Eq.(|8|. Our 

result, however, indicates that as soon as we add in- 
terparticle interactions, the non-linear conductance 
scales as G (X V^~'^ in agreement with the scaling 
theory. Moreover, one can notice that the non-linear 
conductance (26) changes sign across K = 1 which 



is associated with the change of the local density of 
states dependence Pldos(^ = 0,cc;) cx uj^^^~^ [40] ^ 



The results for the tunneling conductance (25) 



and (26) are valid as long as T > A* and eV > A*, 



respectively. Here the scale beyond which the per- 
turbation theory breaks down is defined as A* = 
A(Am/A)^/'^^~^ \ Thus, in order to access the low 
temperature/voltage behavior of the conductance. 
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TABLE 1. Tunneling conductance near perfect normal reflection fixed point, see Eqs 25 and (26) for details, 
we need to compute transport starting from the low- energy perfect Andreev reflection fixed point. 



B. Conductance near prefect Andreev reflection 
fixed point. 

In this section, we discuss transport properties of 
the TP SC/LL junction in the low-energy limit, i.e. 
near the perfect Andreev reflection fixed point. For 
pedagogical reasons, we consider here a simplified 
setup where L ^ oo and relegate the discussion of 
the leads to the next section. We assume that the 
system renormalizes to the perfect Andreev reflection 
fixed point (i.e. the field is pinned), and the effec- 
tive low-energy theory is determined by the following 
imaginary-time boundary action [40]: 



S = 



(28) 



The second term here represents the leading irrele- 
vant perturbation around this fixed point which pro- 
motes normal reflection at the junction. 

The corresponding partition function for the 
boundary theory on a Keldysh contour reads (see 
Appendix [A| for details) 



X exp 



-So[^el,^q]+iSv[^cl,^q] 



(29) 



where [$ci , ^q] is defined as 

So[<^cl,<^q]= I ^($d,*q)c.GilH(<&el,*q)-c. 

(30) 

with G^ ^{uj) being the 2x2 matrix defined in 
Eq. ( |A5[ ). The action Sv[^ch^q] describes coupling 
to the electric field 



Sv[^ci,^q] = -7r / dt[^g{t)Vd{t) + ^ci{t)Vg{t)] 

(31) 



with V{t) being the applied voltage. The last term 
in the exponent of ( 29 ) represents normal backscat- 



tering at the junction 



-2AAr J dt sin ^ci{t) sin ^q{t). (32) 



Given that irrelevant coupling Aat is small, one can 
compute the response function function perturba- 
tively. 

The dc tunneling conductance G can be obtained 
by first computing the current through the junction 
using Eq. (19) and then differentiating it with re- 



spect to the dc voltage V 
d{Ici(t)) 



dV 



^e• 
2^ 



dt'{dt^eimq{t')) 



K-i-0 



(33) 



where the average is defined with respect to + 
Sv^Sxn ' We now perform a transformation of vari- 
ables to eliminate the term linear in V in the action 
by introducing a new variable ^ci{t) = *^ci(^) + <^{t) 
with 



dta{t) = -eKV. 



(34) 



The correlation function in Eq.(33) can be computed 
by doing perturbation theory in A at: 



dt {^el mq {t'))s^,^^,^dt {^el (t) % it')) So (35) 

1 

-2' 



;dt{{^eim,{t')Sl)so- 



Here we used the fact that {a{t)^q{t')) Sq = 0. The 
zeroth order in the Aat correlation function can be 
easily computed using Eq.( |B4| ) yielding the following 
contribution to the dc conductance 



h 



(36) 



In the non-interacting limit K = 1^ we recover the 
previous zero-temperature results ([9|. However, in 
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general Eq. ( 36 ) suggests that conductance depends 



on the Luttinger parameter K. This is a rather sub- 
tle issue [49H54] since the result for G depends on the 
experimental setup. As we show below, the presence 
of the leads is important for the transport measure- 
ments [[UlSHS] and the conductance G = 2e^/h re- 
gardless of the interactions in the LL. 

We now compute the second-order in A at correc- 
tion to the conductance and define the following cor- 
relation function C^{t — t') 

C^{t-t') = \dt{^,im,{t')Sl)s,. (37) 

After some manipulations, one finds that the contri- 
bution of the connected diagrams to the correlation 
function C^{t — t') is given by 



X {5^[-^^q{t'))cos[a{ti)-a{t2)]sm{K2[ti-t2]). 

(38) 

Here = $(ti) ±$(^2); the functions W{t) and 

K2{t) are defined as 



i/t2(ti-t2) = (<5$ir^<5$(+)). 



(39) 
(40) 



The analytical expressions for these functions are 
given in Appendix [BJ After straightforward manip- 
ulations, we arrive at 



C^{t-t') 



-iX 



I 



duj 
2^ 



uit-t'). 



where the function F{uj) reads 

/oo 
d{ti-t2)sm[oj{ti-t2)]e- 
-00 

X cos[a(ti) - ^(^2)] sin[/^2(^i - ^2)]. 



(41) 



Witi-t2) 



(42) 



According to Eq.([34|), the function a{ t) d epend s lin - 
early on t in the dc limit. Using Eqs.(B5) and (B6), 
the function F{uj) is simplified to 



7rT\ 



AK 



dx- 



sinh^^(x) 



X sin(27ri^) cos 



eVK 



(43) 



One can notice that the integral (43) is UV-divergent 



in the regime of interest, i.e. at > 1/2. This is 



a common issue, often appearing in boundary sine- 
Gordon quantum impurity problems [68^, '69 when 
one considers fiowing back from the IR fixed point. 
One way to proceed here is to employ dimensional 
regularization [69|, i.e. we compute the integral in 
the domain of where it is well-defined, and then 
analytically continue to K > 1/2. 

At small frequencies uj <^T^ one can simplify the 
expression for F{uj) 



F{u;) 



2ui /'2iT\4K-2 
A2 I A j 

A2 I A j 



gi{K), eV/T^l 
g2{K), eV/T » 1. 



(44) 



Here the dimensionless functions gi{K) and g2{K) 
are given by 



gi{K) = -2^^-^itcos(2ttK) 



T{-4.K)T(l^2K) 



g2{K) = -2K{4K 



r(l - 2K) 

(45) 

l)r(-4i^)sin(47ri^). (46) 



In the non-interacting limit K = 1, these functions 
are finite: gi{K = 1) « 7r/3 and g2(Jf = 1) ^ 7r/4. 



Taking into account Eqs.(33), (41) and (43), the 



corrections to the dc conductance can be now calcu- 
lated in both limits: 



SG^'^ = 



h 7rA2 



tA2 



{^y—'g2{K) eV»T. 

(47) 



As one can see, the scaling of the dc conductance is 
consistent with the non-interacting results discussed 
in Secini 



C. Conductance near prefect Andreev reflection 
fixed point: eff'ect of the leads 

In this section, we consider a finite-length LL and 
discuss the effect of the leads on the conductance 
near the low-energy perfect Andreev reflection fixed 
point. The experimental setup we consider here is 
shown in Fig{l] We will use here an inhomoge- 
neous Luttinger-liquid model and introduce an in- 
homogeneous interaction parameter K{x) and ve- 
locity v{x)^ which adiabatically interpolate between 
the nanowire and lead, see Ref. for more de- 
tails. Since electrons in the lead form a Fermi liquid, 
which for the present purposes is treated as a non- 
interacting Fermi gas, we assume here that Kl = 1. 
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In this case, the naive RG analysis does not work 
due to the spatial dependence of K{x) and v{x). 

We, therefore, pursue a different route here and 
first compute the Green's function for the LL liquid 
coupled to the leads at x = L/2 and the TP SC 
at X = 0. At the perfect Andreev reflection fixed 
point, 0{x = 0,t) is pinned, and, thus, an appro- 
priate boundary condition is dtO{x = 0, t) = 0. In 
terms of the fluctuating variable ^, this is equiva- 
lent to dx(l){x = 0,t) = 0. Similarly, one can define 
boundary conditions at x = see Appendix [Cj 
The real-space Green's function for the Luttinger 
liquid coupled to the leads can be found by solving 



UJ 



^ K{x)v{x) 



)} 



Guj{x,x) = 5{x - x) 
(48) 



where we take the simplest ansatz for the spatially- 
inhomogeneous velocity v{x) and Luttinger param- 
eter K{x) [49] 



v{x) 



K{x) 



V x<L/2 
v' x> L/2 

Kw x< L/2 
Kl x>L/2. 



(49) 



(50) 



The straightforward solution of the above wave 
equation Goj{x,x') at the origin x = x' = {) (see 
Appendix [C| for details) reads 

One can notice that there is an additional energy 
scale in the problem, = v/L, which is associated 
with time of propagation through the nanowire to 
the leads. 



Using Eq.(51), it is straightforward to obtain the 
modified Keldysh action in the presence of leads at 
the Andreev fixed point. 

5o[$d,*,]= j ^(*cl,*q)wGil(w)(*cl,*q)-o. 

(52) 

where the Green's functions are given by 



'oj±i5' K^+KL-{KL-K^)e±i-^L/v 

(53) 

27TiK,„ / a; \ 



coth I - 



\2TJ 
2KlK^ 



(54) 



Kl^Kl-{Kl-Kl)co^{^y 



We now compute the conductance in a similar way 
in Sec. [mB) The full Keldysh action is Sc = 
^Xj^ where So^Sy and are defined in 



as 

Eqs.(l52l), (|3l]) and (|32]), respectively. Using Eq.Q, 



we first compute the conductance to zeroth order 
in A AT. It is straightforward to show that the dc 
conductance G^^^ (cf. 36) is given by 



i—r-^GS{^)\ou^O 
TTfl 



2e^ 



Kr 



(55) 



Note that the conductance depends on Kl = 1 
rather than K^, which is consistent with the pre- 
vious works on a similar problem j49H5T| [53l 154] , 
Thus, in a typical setup for measuring dc transport, 
see FigjT] the perfect Andreev reflection fixed point 
corresponds to a quantized conductance in units of 
2e^ /h even in the presence of the interactions in the 
nanowire. 

We discuss next the corrections to conductance at 
finite temperature and voltage. In order to compute 
the corrections to the lowest non- vanishing order in 
Aat, we first shift the field ^ci{^) <I>cz(t) +a(t) and 
eliminate the Sy term in the Keldysh action. In the 
case of an applied dc voltage V, the function a{t) is 
given by 



a{t) = 2eV 



duj 
2^ 



-iujt 



Using Eq.(53), one finds that dta{t) = -cVKl. The 



rest of the calculation proceeds exactly as in the ho- 
mogeneous LL case discussed in Sec |III B The only 
difference here is that the correlation function (37) 



is now averaged with respect to the action 5*0 (52): 



C^(t-t') =Ai\% 



I 



duo 



-iuj{t-t') 



(56) 



The function F{uj) is given by 

/CO 
dr sin(cL;r)e~^'^'^^ cos[eVLT] sin[/^2(T)] 
-oo 

(57) 

where W{t) and /^2('?") are defined in Eqs.(|D2[) and 



(D3), respectively. In general, W{t) and h2{,^) are 
complicated functions of Ky^ and Kl^ and the in- 
tegral (57) has to be evaluated numerically. One 



can also evaluate these functions analytica lly i n the 
limit of small 5K = K^^ — Kl, see Eqs.(D2) and 
(|D3l) 



which allows one to compute the conductance 
in the limit of weak electron-electron interactions in 
the nanowire. Henceforth, we assume that \5K\ <C 1. 
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The function F{uj) also depends on the ratio be- 
tween T, eF, and El, i.e. the main contribution to 
the integral ( [57| ) can come from the interval r < L/v 
or r > L/v subject to the infrared (IR) cutoff de- 
termined by max(ey,T). If the wire length is much 
larger than the thermal coherence length Lt, the 
correlation function F{uj) is approximately given by 



F{uj) 



2 /ttT 



dx 



X sm{2TTKw) cos 





cVKl 



sinh^^-(a;) 



(58) 



At small frequencies uj <^T^ one finds that 



F{u;) 



A2 



(if) 



4K^ 



2u; [ eVKL 
A2 V A 



'gi{K^), eV/T^l 

~%2{K^), eF/T»l 

(59) 



where the functions gi{K) and g2{K) are defined in 
Eqs.(|45|). 



In the opposite limit L <^ Lt, the function F{uj) 
is approximately given by 



/CO 
:T1 



2 



sinh^^^(x) 



(la; 

sin(27r i^T/,) cos 



(60) 



cVKl 
ttT • 



In the low-frequency limit co -t^T, the expression for 
F{lo) simplifies to 



F{oj) 



2uj ( V 
A2 VlA 



A5K 



[ ttT 
\ A 



200 f V 



(fx)"" (^) 



A^ 



4Kr 



'5i(ifL), 



■52(^^1 



eV 



> 1. 

(61) 



Using Eqs.(33),(56) and (l57|), we can readily ob- 



tain the dc conductance in the presence of leads. As 
expected, in case oiL^ Lt the presence of the leads 
is irrelevant, and we obtain the same answer as in 
the L ^ 00 limit, cf. Eq. (47) where K is the Lut- 
tinger parameter in the wire. On the other hand, if 
L <^ Lt one finds that the temperature dependence 
of the corrections to G is governed by the Luttinger 
parameter of the leads Kl. The conductance G in 
this regime is given by 



G 



G^ 



2e' 



h A2 V A y 



2e2 



-Kl- 



8X%K 
" h A2 



-2) 



(^) 



4SK 



gi{KL) for eV <.T 



(62) 



A ; 



-2) 



Ok) 



4SK 



g2{KL) for eV>T 



r 



where we eventually have to set Kl = 1. Thus, 
in this limit temperature/voltage corrections to the 
tunneling conductance scale similarly as in the non- 
interacting case. One can notice, however, that SG 
has a non-trivial dependence on the nanowire length 



L and the Luttinger parameter K^. In the case of 
repulsive interactions in the LL {K^ < 1), the cor- 
rection to the conductance is enhanced by the factor 

Kla) 



r 



IV. QUALITATIVE DISCUSSIONS OF 
MAIN RESULTS 



ergies (i.e. near the perfect normal reflection fixed 



We now review our main results for the tunneling 10 
conductance in the TP SC/LL junction. At high en- 



eV <^T 



e2 8Aj^ ^^T\(^K^-2) , 



6G^ 



2) 



L<^Lt 



"TT l^J [la) 9i[J^l) 



eV^T 



2 A j 92[JAu 



h A2 



^ A2 



^1 f eVKL \ 



TABLE II. Corrections to the tunneling conductance near perfect Andreev reflection fixed point, see Eq.(47) and (62) 
for more details. Here Lt = v/T is the thermal coherence length. 



point), the results for the tunneling conductance G 
are in Table 1. We find that the temperature and 
voltage dependence of the conductance are consis- 
tent with scaling theory. However, the pre-factor of 
the non-linear conductance ( 26 ) vanishes in the non- 
interacting limit = 1. This is consistent with 
the previous results based on the non-interacting 
calculation, see Sec. Thus, the non-interacting 
result for conductance is quite fine-tuned and does 
not represent the generic situation in the nanowires 
where interactions are inevitably present. We show 
that the non-linear conductance changes sign across 
the = 1 point which is associated with the 
interaction-induced enhancement /suppression of the 
local density of states at the junction. 

The low-energy results for the tunneling conduc- 
tance G obtained near the perfect Andreev reflec- 
tion fixed point are summarized in Table II. First 
of all, we showed that in the presence of the leads, 
the dc tunneling conductance is equal to 2e^ /h at 
zero temperature. This is consisten t with the pre- 
vious studies [49H5n [53] . In Sec. |IIIB| we stud- 



ied finite-temperature and finite- voltage corrections 
to the tunneling conductance and showed that the 
power-law dependence on T and V is not universal 
and non-trivially depends on the length of the LL re- 
gion L, see Table II. These results can be understood 
from a very general perspective using a two-step RG 
analysis. 

Let us consider the scattering from the bound- 
ary at X = 0. In the unfolded picture discussed in 
Sec. |TI| the boundary can be viewed as an effective 
"impurity". It is well known that charge density 
oscillations (i.e. Friedel oscillations) around such 
an "impurity" lead to the dynamically-generated 
effective potential which determines the transmis- 
sion/reflection probability. The spatial extent of the 
Friedel oscillations is determined by the thermal co- 
herence length Lt = v/T. Thus, at high tempera- 
ture Lt <C the reflection probability from the im- 
purity depends on v and i^^, whereas in the opposite 



limit the effective size of impurity potential becomes 
larger than the length of the wire, and scattering is 
governed by the dynamics in the leads characterized 
by v' and Kl- These arguments suggest that one 
can understand the results (62) using a two-step RG 



flow of the coupling constant A at. In the homoge- 
neous case (i.e. L oo), the RG equation for A at is 
given in Eq. ( 15 ) where K is the LL parameter in the 
nanowire. If the LL length L is finite, there is an ad- 
ditional energy scale in the problem El = v/L which 
separates the two aforementioned regimes. Specifi- 
cally, in the high-energy (short-length scale) regime 
A > £^ > El^ the flow of the coupling constant Xn 
is governed by K^j. If max{T, eF} El, we have 
to stop the RG flow before it reaches El- Therefore, 
the expression for the effective coupling constant at 
the IR cutoff scale is independent of L: 



2K^ 



The correction to the conductance is proportional to 
(A^)^ and is given by 



JG(x-A^(0) 



2 . _ . / max{T, eV} 



A 



(63) 



In the case max{T, eF} <C El, the renormalization 
of Xn is determined by and Kl in the energy 
intervals A > E > El and El > E > max{eF,T}, 
respectively. A straightforward integration of the 
RG flow for A AT yields 



= A^(0) 



2K^-l 



max 



{eV.T} \ 
El ; 



and one finds that the correction to the conductance 
is given by 



(64) 



11 



The scaling of conductance ( 63 ) and ( 64 ) agrees per 



fectly with the results obtained using Keldysh for- 
malism, see Table 11. 

Thus, we have shown that the two-step RG ap- 
proach significantly simplifies the calculation of tem- 
perature/voltage corrections to the tunneling con- 
ductance and allows one to take into account the 
effect of the leads. This approach is quite general 
and can be easily applied to other quantum impurity 
problems where the presence of leads is important. 
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Appendix A: Conventions and Notations 



V. CONCLUSIONS 



In this paper, we considered transport proper- 
ties of a topological superconductor-Luttinger liquid 
junction and computed the tunneling transport con- 
ductance using a real-time Keldysh technique. We 
have considered a realistic experimental setup [l] 
and taken into account the effect of the leads within 
an inhomogeneous Luttinger liquid model. We find 
that interactions in the Luttinger liquid do not mod- 
ify the quantization of the zero-bias tunneling con- 
ductance at zero temperature but rather lead to a 
change of the temperature/ voltage corrections to the 
conductance. We show that in some instances the 
results based on non-interacting models are signif- 
icantly modified due to the presence of the inter- 
actions. Thus, our results have important impli- 
cations for the tunneling experiments aimed at de- 
tecting Majorana zero-energy modes in semiconduc- 
tor / superconductor heterostructures. 



In this Appendix, we specify our conventions and 
review basic properties of the Keldysh technique 
used throughout the main text. We follow closely 
the Keldysh path integral approach reviewed in 
Ref.[67 . We start with the definition of the par- 
tition function on a Keldysh contour which reads 

Z = J D(l)exp(^i J dt(t)(t)G-^(t)(t)^ . (Al) 

The fields on the upper and lower branches of the 
Keldysh contour are ^+ and respectively. In 
general, the action Sc can be written as 

Sc= I dt(t){t)G-^(t){t) = [ dt(t)^{t)G-^(t)+{t) 

Jc J-oo 

oo 



Next, we perform a Keldysh rotation by introducing 
new fields 0ci = ^++^- and 0q = — Without 
loss of generality, we now consider the action for the 
finite-length Luttinger liquid, which, after Keldysh 
rotation, reads 



rL/2 

Sc= I dt I dx 



o(nci,nq) 



dt 
dt 



1 / \ V 



(A2) ■ 



r 



By integrating out lici- and Ucrfields in Eq.(A2), where the matrix Green's function ^ is defined as 
one finds 



= J D(t)c\D(t)^ exp 



(A3) 



x,t J x' ,t 



(A4) 



G-/ 



(A5) 
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The Green's functions are given by 



5{x - x')8{t - t') (A6) 



2t:Kv 

and {G~^)k is defined as [67] 

{G-^)K = GR^oF-Fog-\ (A7) 
In the translationally-invariant system we find 

/ N 27TVK . , , 

g^/^ q) = ^ A8 

q) = coth (^) q) - g^ioj, q)] 

' uj \ ttK 



= coth(^- 



—2m ^{uo^vq)-5{uo- vq)] . ( A9) 



The following correlation functions are defined in 
terms of the Green's functions as 

{^^i{x,t)^q{x' ,t')) = 2ig^{xt',x't') (AlO) 
{^q{x,t)(j)ci{x',t')) = 2ig^{xt',x't') (All) 
{Mx,t)(})ci{x',t')) = 2ig''{xt;x't'). (A12) 

We can now derive the bound ary t heory by inte- 
grating out the bulk modes in Eq.( A3 ). The Keldysh 
action with the constraint ^{t) = (j){x = 0, t) is given 

by 



D<pciD(f>gDxciDxq exp 



+ 1 J dtdx {xq{t)5{x)[^ei{t)-(i)ci{x,t)] 
+Xci{mx)[^S)-<P<i{^M\- (A13) 

We first integrate out the (l){x^t) fields and then per- 
form a functional integral over xifi)- As a result, one 
arrives at the following boundary theory 



ZcB = j D^^iD^qexp -Sol^ci.^q] (AM) 
So[^ci.^q]= (A15) 
{^,u%)tG-\x = 0,t;x' = 0,tO(^ci,^q)^'. 



This result is quite general and applies to the inho- 
mogeneous Luttinger liquid model discussed in the 
main text. 

The correlation functions for the ^-field can be 
straightforwardly obtained by changing K 
in Eqs.(|A6|) and (IA7). 



Appendix B: Correlation functions for 
semi- infinite LL 



In this section, we derive the expressions for the 
various correlations functions for the semi-infinite 
LL case. Due to the translational invariance in the 
wire, the Green's functions depend on the difference 
x — x'. Therefore, the retarded and advanced bound- 
ary Green's functions become 



I 



duj dq 



27T 27T {vq — uj — i5){vq -\- uj -\- i5) 

(Bl) 



where a is a cutoff. In the frequency domain, we 
find 



iriK 
{uj ± iS) 
2tiK 



coth {^) 



(B2) 
(B3) 



with corresponding real-time Green's functions 



irK 



27r u ±i5 

= -iTKd{±M) (B4) 

where ^{t) is the Heaviside step function. We note 
that the equal-time Green's function is equal to zero 
5'o[0c^0g] (0) ~ 0' seen by performing the fre- 



quency integral first in Eq.(Bl). 



Let us now compute several correlation functions 
that appear in the calculations. The function W{t) 
appears in the perturbative expansion 



W{ti-t2) = 

POO 

= 4K dx 
Jo 



(($e/(il)-*c/(i2))') 
- COS 



= 8K 



dx I 1 — cos I X 



ttT 



coth (I) 



exp[ 



A^ 

= 2i^log( -^sinh'[^T(|ti-t2| 



(27rn)2+x2 
(B5) 



where a = A/T with A being an UV cutoff. 
Also, we have assumed here that \ti — ^ A~^. 
As can be easily shown, the correlation function 
{{^ci{h) + *^d(^2))^) is divergent and, thus, the 
terms in the perturbative expansion proportional to 
exp[-^{{^ci{ti) + ^d(t2))^)] vanish. 
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The correlation functions Ki{ti—t2) and ^2(^1—^2) 
are defined as 

= 2TTiK {^{h - t2) + ^{t2 - h)) 

= 2mK (B6) 

iK2{ti - 12) = (.5$^-^,5$(+)) = -27riii:sgn(ti - ta) 

(B7) 

where 5¥^^ = <^>{ti) ± ^(ta)- 

The correlation functions for the ^-field can be 
similarly obtained. Using the duality between the 9 
and (h variables, one finds 



^imir) = ^log ( ^sinh^[7rT(|ti-t2| 



{sei;^seW) = -^-^sgn{t,-t2). 



(B8) 
(B9) 

(BIO) 



Appendix C: Green's function for LL coupled to 
the lead. 



We now consider the effect of the leads and derive 
the real-time Green's function for a LL of length 
L/2 coupled to the metallic lead. We will use here 
the inhomogeneous LL model [49] and first solve the 
following wave equation for the Green's function in 
the frequency-position domain 



_ K{x)v{x) 



(Cl) 



We take the inhomogeneous velocity v{x) and Lut- 
tinger parameter K{x) to be 



v{x) 
K{x) 



V X <L/2 
v' x> L/2 

X < L/2 
Kl x>L/2. 



(C2) 
(C3) 



Since K{x) and v{x) are constants at each point in 
space, the Green's function takes the general form 

= A(xOe^^^ + 5(xOe"^^^. (C4) 
The Green's function should satisfy continuity at 
X = L/2 as well as the following boundary condi- 
tions at X = 0, x', and L/2: 



v{x) 



K{x 
v{x) 



dxGuj{x,x) 



x=0 



K{x) 
v{x) 



dxQuj{x,x') 



K{x] 



dxQuj{x,x') 



=L/2 
x=x' 






= 1. 



(C5) 
(C6) 

(C7) 



We are ultimately looking for the solution for the 
Green's function in the wire and henceforth assume 
x' < L/2. This divides our space into 3 regions, 
where we must match the 4 boundary conditions 
above 

Q^^\x,x') = Ai(xOe^^^ + Bi{x')e-'^^ < x < x' 
gS^\x, x') = A2{x')e'^'' + B2{x')e-'^'' x' < x < L/2 
= A:,{x')e'^^ X > L/2. 

Solving for the coefficients Ai and Bi, we find the 
Green's function at x < L/2. 



iK^ eiL/L^^2_^-i\x-x'\/L^ ^ ^-iL/L^^2^^i\x-x'\/L^ ^ 2/^_/^+ COs(|x + x'\/L^) 



where i<i± = \/Kyj ± and L^^ = v/u. Upon 

setting X = x' = 0^ one arrives at Eq.(53). 



Appendix D: Boundary correlations Green's 
functions for the inhomogeneous Luttinger 
liquid model. 



We now compute the relevant correlation func- 
tions in the presence of the leads. In general, one 
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has to compute these numerically. One can obtain, 
however, some limiting expressions for \5K\ <C 1, 
where 5K = — K^. We first compute the corre- 



lation function W{ti — ^2) 

W{h-t2) = \{[^ci{h)-^,i{h)f) (Dl) 

/OO 7 
^{l-COs{L0[h-t2]))2ig^{L0). 
-00 

By expanding (oj) to first order in SK, the corre- 
lation function W{At) becomes 



W{h-t2) 



— (1 — cosiujAt)) 



+ cos 



2T) 



: 2K^Y{M) - SK [Y{At + L/v) + Y{At - L/v) - 2Y{L/v)] 

I 



(D2) 



where Y{At) = log (^;^sinh^[7rTAt]y Here we 
implicitly assumed that \At\ A~^. 



The correlation functions 1^1 and 1^2 in the pres- 
ence of the leads are given by 

ik,{At) = {S^[i^S^[-^)=4i rp[l-cos{ujAt)]g§{uj) 

5K 
(D3) 

iK2(Ai) = (5$(-)5$(+))=4 r ^sm{u;At)g§{uj) 

J -00 









\^ . 1 


^2 — sgn 




— sgn 













= — 27rzi^r^sgn(At)— Tri^i^l sgn 









X , " 




^sgn 




— sgn 




) 




V 




V 





(D4) 
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